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Introduction. 

Clebsch has shown that the components of the velocity of a fluid u, v, w, parallel to 
rectangular axes %, y } z, may always be expressed thus 

U — -f-A — , V — ~tA~, W=—-t-A ~^r I 

ax dx ay ay dz dz 

where X, \fj are systems of surfaces whose intersections determine the vortex lines; and 
the pressure satisfies an equation which is* equivalent to the following 






dx) \dy) \dz 



+w 



\dx ) \dy } \dz ) 



where p is the pressure, p the density, and V the potential of the forces acting on the 
liquid. 

It is shown in this paper that an equation of a complicated nature in X only can be 
obtained in the following cases (that is to say, as in cases of irrotational motion, the 
determination of the motion depends on the solution of a single equation only) ; — 

(1.) Plane motion, referred to rectangular coordinates x, y. 

The equation is somewhat simpler when the vortex surfaces are of invariable form, 
and move parallel to one of the axes of coordinates with arbitrary velocity. 

(2.) Plane motion, referred to polar coordinates r, #. 

The equation is somewhat simpler when the vortex surfaces are of invariable form, 
and rotate about the origin with arbitral angular velocity. 

(3.) Motion symmetrical with regard to the axis of z in planes passing through it, 
referred to cylindric coordinates r, z. 

The equation is somewhat simpler, when, the vortex surfaces are of invariable form, 
and move parallel to the axis of z with arbitrary velocity. 

* British Association Report for 1881, p. 62, 
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Suppose that in any of these cases any particular integral of the equation in X is 
taken. 

It is shown that the components of the velocity can be expressed in terms of X and 
differential coefficients of X, and that the current function is also known. 

In the case of a fluid, part of which is moving rotationally and part irrotationally, 
the boundary surface separating the rotationally moving fluid from that which is 
moving irrotationally contains the same vortex lines, and may be taken at the surface 
X=0. 

Now, if the integral taken of the equation in X do actually correspond to a case of 
fluid motion in which part of the fluid is moving rotationally and part irrotationally, 
the most obvious way to find the irrotational motion will be to find its current function 
from the conditions supplied by the fact that the components of the velocity are con- 
tinuous at the surface X=0. Examples I. and III. of this paper have been solved in 
this manner. 

If after taking any integral of the equation in X it be found theoretically impossible 
to determine the current function of an irrotational motion outside the surface X=0, 
which shall be continuous with the rotational motion inside it, then the integral in 
question does not correspond to such a case of fluid motion. 

In this method no assumption is made as to the distribution of the vortex lines (as 
in the method of Helmholtz) before commencing the determination of the irrotational 
motion. 

If, however, the rotational motion be known, the components of the velocity are 
known for this part of the fluid. Let the components of the velocity be expressed in 
Clebsoh's forms, so that x> X, \jj are known. 

Moreover, let the forms be so arranged that the surface separating the rotationally 
moving fluid from that which is moving irrotationally is the surface X=0. 

^i i •. dy dy dy 

Then at this surface the components of the velocity are -~> -~> -£• 

1 J dx dy dz 

Now, obtain in any manner a velocity potential <f) for space outside X=0 continuous 

with x a H over the surface X=0. This is theoretically possible always. 

If the velocity potential so obtained make the velocity and pressure continuous all 

over the surface X=0, then a possible case of motion will have been obtained. 

The conditions to be satisfied in order that the velocity may be continuous at the 

surface X=0 are that there ~^=-^, J^=^, ^= J_. In order that the pressure may 

ax dx dy dy dz dz 

be also continuous, it is further necessary that ~~=~^ all over the surface X=0. 

The most obvious way of obtaining the velocity potential will be to apply Helm- 
holtz's method of finding the components of the velocity in terms of the supposed 
distribution of magnetic matter throughout the space occupied by the rotationally 
moving fluid. 

It must, however, be remembered, as is remarked by Mr. Hicks in his report to the 
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British Association on " Recent Progress in Hydrodynamics/', Part 1,* " That the 
results refer to the cyclic motion of the fluid as determined by the supposed distribution 
of magnetic matter, and do not give the most general motion possible/' It appears 
also from Examples I. and III. of this paper that it is not possible to assume arbitrarily 
the distribution of vortex lines, even when it can be shown that the equations of 
motion are satisfied at all points where the fluid is moving rotationally, and then to 
proceed to calculate the irrotational motion by means of the supposed distribution of 
magnetic matter. For in these examples, values of the components of the velocity of 
a rotational motion, satisfying the equations of motion throughout a finite portion of 
the plane of x, y, are found. Thus the distribution of vortex lines, and, therefore, that 
of the supposed magnetic matter over a finite portion of the plane of x 9 y is known. 
The surfaces that always contain the same vortex filaments are found. Inside one of 
these the supposed magnetic matter is distributed, the current function at external 
points is calculated by Helmholtz's method, and it is shown that the velocity thence 
deduced is not continuous with the velocity of the rotational motion at the surface, 
which separates the rotationally moving liquid from that moving irrotationally. 

Another way (suggested by Clebsch's forms) of obtaining the velocity potential will 
be as follows : — 

Calculate the quantity p=~-t~(-t^+^+-t? 

Treating p as the density of a material distribution inside k=0, taking no account 
of the value of p outside the surface X~0, obtain the potential of this distribution. 
Let the potential inside X=0 be x'> and outside let it be <£. 

X will, in general, differ from x 5 fast, because x ma y contain many-valued terms, 
which may be denoted by P, satisfying Laplace's equation ; and, secondly, because 
X— P may be the potential of a distribution of matter, part of which is outside X=0. 

Accordingly, it is necessary to examine whether it is possible to find many-valued 
terms P satisfying Laplace's equation such that x'+P^x- 

Then (£+P will be the velocity potential of the irrotational motion, provided that it 
give zero velocity at infinity. 

Example II. of this paper is solved in this manner. It might also have been solved 
by Helmholtz's method. 

The few illustrations which follow are a first attempt to apply the theory to 
particular cases. 

Example I. treats of the motion of an elliptic vortex cylinder of invariable form 
parallel to one of its axes with arbitrary velocity. The irrotational motion outside 
the cylinder cannot be supposed to extend to an infinite distance. 

Example II. treats of Ktrchhoef's elliptic vortex cylinder, in which the angular 
velocity of the rotation of the cylinder is a function of the vortex strength, and the 
axes of the elliptic section of the cylinder. 



* 
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Example III. treats of the revolution of an elliptic vortex cylinder round its axis, 
where the angular velocity is not restricted as in the last case. The irrotational 
motion outside may be supposed to be limited by a smooth rigid confocal elliptic 
cylindric surface, rotating with the same angular velocity. The last example is the 
particular case of this, obtained by supposing the elliptic section of the external 
confocal cylinder to become infinite. 

Example IV. treats of the motion of the fluid in a fixed circular cylindric surface, 
where the vortex strength is any function of the distance from the axis, the irrotational 
motion continuous therewith being supposed to extend to an infinite distance. 

Example V. treats of a possible case of rotational motion inside a certain hollow 
smooth rigid surface of annular form, which moves parallel to its straight axis with 
arbitrary velocity. 



1. Clebsch's * forms for the components of the velocity of a liquid u, v, w parallel to 
fixed rectangular axes x, y y z in space are : — 

%l — — -f-A— , V—- \-k-r~ , IV — — -f-A , 

ax ax ay ay dz az 

where the surfaces X= const., i/f= const, determine by their intersections the vortex 
lines, and always contain the same particles of liquid. 
If F(X, \p) be an arbitrary function of X, xjj 

ad/ adj 

and similar expressions for v, w. 

Thus these expressions for the components of the velocity are still in Clebsch's 
form. 

\\ x\j are each functions of X, xjj. 

X satisfies the same equation as x- 

* Taking as independent variables three families of surfaces, always containing the same particles, and 
the time, the writer obtained independently Clebsch's forms in an article published in the Quarterly 
Journal of Pure and Applied Mathematics, February, 1880, vol. xvii., entitled " On Some Properties of 
the Equations of Hydrodynamics." 

A demonstration of the same forms for any fluid in which the density is any function of the pressure 
is contained as a particular case in a paper entitled " On Some General Equations which include the 
Equations of Hydrodynamics," which is published in the Transactions of the Cambridge Philosophical 
Society, vol. xiv., part i., the writer having previously seen Clebsch's paper, " Ueber Die Integration der 
hydrodynamischen Gleichungen," ' Crelle,' Bd. lvi., p. 1. 
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Moreover, since F(X, xjj) is an arbitrary function of X, xjj ; it can be so chosen that X' 
may be any required function of X, \ji ; i.e., any vortex sheet. 

Therefore the X in Clebsch's expressions for u, v 9 w may be considered as the 
surface of any vortex sheet ; and^ consequently, as the surface separating the 
rotationally moving fluid from that which is moving irrotationally. 

Therefore, if <j> be the velocity potential of the irrotationally moving fluid, all over 
the surface X=0, supposing the motion continuous there ; 

d% d<f> dy dcj> dy d<f> 

dso dot ' dy dy 9 dz dz' 



Since also 






n + V = 
P 


d X 

dt 


i (f d x 

2 \\dx 



'dfy p fdf\* fdfY 

V dz 



in the rotationally moving fluid, and 



p dt 2 \\dxJ ' \dy) l \dz / 

in the irrotationally moving fluid, it follows that the condition for the continuity of 
the pressure at the surface separating the rotationally moving fluid from that moving 

irrotationally is that -~ = -^ all over the surface X=0. 

Now suppose that there exists a solution of the three equations to w^hich x, X, \p are 
subject as given by Clebsch ; then to find <£, it is necessary to find it so as to satisfy 
the above surface conditions. 

In any case in which x is the potential of a distribution of matter inside the surface 
X=0, together with many valued terms satisfying Laplace's equation, then <f> is the 
potential of this distribution calculated for a point outside X=0, together with the 
same many valued terms, provided that it give zero values for the components of the 
velocity at infinity. 

With regard to the supposed distribution of matter, its total mass must be zero, in 
the case of an incompressible fluid. 

For total mass of supposed distribution of matter = — — ~ dS, the integration 

being extended over the surface X=0, 

1 
But this =— — (total flux outwards across the surface) = 0. 

4/r x ; 

Therefore total mass of supposed distribution of matter =0. 
2. Plane Motion. Rectangular Coordinates. 
In the ordinary notation the equations are 
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du 
dt 

dv 
dt 



dib du 



dx 



dp 



-p V 



dv dv d /fdp ~ 

dx dy dy\) p 



j;+tAp u )+^ / (p v )=° 



dt dx 



From which can be deduced 



d 



d 



d 



dt ' dx dy 



< 



dv du 
dx dy 



> = 



v». 



Now regarding this as a particular case of motion in three dimensions in which 
tvz=0, the motion being parallel to the plane z=0, it is possible to put by means of 
Clebsoh's forms 



Therefore 



u== d2c x a± m v== dx d± 

dx dx ' dy dy 



1 J J \ f^d^fr dX difr~) 

-f^~K dx dy dy dx > = 



dt dx~ 1 ~ l dyj 



This result has been deduced from the three equations of motion, and Clebsch's 
forms for the components of the velocity. 

But it can be deduced from the equation of continuity alone and the following 
equations known to be satisfied by X, %jj. 



Therefore 



u= 



dX dX dX 

dt dx dy 







dty d^Jr ^__ n 

dt dx dy 



dX 


dX 


dt 


dy 


dty 


dyjr 


dt 


dy 


dX 


dX 


dx 


dy 


dty 


dty 


dx 


dy 



and v- 



dX 

dx 


dX 
dt 


dyjr 

dx 


dty 
dt 


dX 


dX 


dx 


dy 


d^fr 

dx 


dty 

dy 
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JLnereiore 



dib dv 
dx dy 



dPX_ 
dxdt 

d 2 ty 

CtfXC&t/ 



dX 

dy 

dy}r 

dy 



■MmZm* 



dX 
dt 

dty 
~dt 



dX dX 



(XX 

d^r 



dy 
dty 



dm dy 



d?X 




dX 


dX 


dxdy 




dt 


dy 


d 2 ^jr 




dyjr 


d^r 


dxdy 


<■ 


dt 


dy 






dX 


dX 2 




CvX 


dy 




dty 


dty 






(XX 


dy 



d_ 

(X'X 



dX dX 
dy 

dty 



d'xjr 

dx 



dy 



But from the equation of continuity 





h d 2 X 


dX 




dX 


#X 




dX 


dX 




dxdy 


dt 
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dx 


dydt 




dx 


dt 




d^ 


dyjr 


dyfr 


d?ty 




d-\}r 


d-^r 


4- 


dxdy 


dt 




dx 


dydt 




dx 


dt 


-|- 


dX dX 




dX 


dX 2 




dx dy 




dx 


dy 




dty dyjr 




dty 


d^jr 






(XX 




dy 






dx 


dy 



A. 

dy 



Therefore 



du dv \ Id d d , 

1. 1 "' ' '' ' -<— J ■ I— ^ ' — j 0} I t\ 

dx dy p \dt dx dy; r 



\ Id d d\ 





d*X , 


dX 




dX 


d*X 




dX 


dX 




dxdt i 
d 2 yfr i 


dy 
dty 


■Jm 


CvX 

dty 


dydt 
d^ 




dt 
dyfr 


dy 

d^r 




dxdt < 


dy 




c&x 


dydt 




dt 


dy 




dX dX 






dX 


dX 2 




dx dy 




(XX 


dy 




dty d^r 




dyfr 


dyjr 






CvdD 
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iy 






dx 


dy 



d 

(X/X 



dX dX 



dx 

dyfr 

ax 



dy 

dty 
dy 



dX 

(XX 

dty 

(X/X 



dX 
dy 
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dy 



dX 


dX 


dx 


dt 
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dx 
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dX 
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Therefore 



(d d d\ 

\dt dx dy) " 



1 Id 
P 



d , d . d 

— „ I , fit -JL.fi) — 

M dx dy / 



Ineretore 



But since also 



dX dX 

dx dy 

dx dy 



dk dty dX d*$r 
d , <# , d\dx dy dy dx ^ , ^ 

T.+^T+^rr — ~ : =0 as before. 

dt dx dy p 



I I \ > /N 

<# a$ <% 



dX 

Q/X 

dty 



dX 
dy 

d'xfr 



dx dy 



a 



nd 



~\-u -j~-\-v — ) \jj=0 



d d . d 

dt dx dy 



Therefore 



dX d^fr dx dty 

^ZJOL — ^lJE. = S ome function of K i|r=/(X, i/0 

r 



3. In what immediately follows p will be supposed constant. Using suffixes to 
denote differential coefficients 

Now let g(\, xjj) be a function of X, \jj such that 



ty(X> f ), 



ineretore 



\ bf 



b^r f(\ f) 






^)=i 



therefore 



M fcfr ^ + iX M M fcfr '* + bX A * 



therefore 









Treating this as a partial differential equation for g, the auxiliary system of 
equations is 

dx dy dg 

Xy X x X 



OF WHICH IS MOVING ROTATIONALLY AND PART IRROTATIONALLY. 371 

Whence each of these 

X x dx +■ X y dy 
— - 

Since, in the differential equation, there are no differential coefficients with regard 
to t y X= const, is one integral of the auxiliary system. 

By means of the equation X= const., y can be expressed as a function of the constant, 
x f and t. As the constant will have to be replaced by X afterwards, it may be said 
that y can be expressed as a function of X, x 9 t ; and when this value of y is substituted 

in — \ y , let the result be denoted by (—K) x * Let differentials of the variables X, x, t 
regarded as independent be denoted by 6X, bx, bt respectively. 

Then g= — - J — ^+ an arbitrary function of X, t. 

1 1 will be convenient to write the arbitrary function in the form - rr~ 
Therefore 



9 



bx , i>F(\, t) 



(ky)x 



K+ bX 



h~F(X f\ 

Therefore g will appear in the form G(X, x 9 i)-\- — jz~* 

The equations of the vortex sheets are functions of \fj, X and therefore of g } X. 
Now 

X^+w\j.+ 17X^=0 



and 



therefore 



9t+ug z +vg y *=0 



bt + bx t+ btbx"^ bx* 



that is 



therefore 



But 



therefore 









bt^ btb\ ^ bx 



bx ~Ky 
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substituting this value of u in the equation 



\ t +u\ x +v\ ¥ =0 



it may be shown that 



v- 



* l bt + btbx 



X 



y 



To determine the current function A, there are the equations 



— — \(— h 
dy \ bt, 



te_ b 8 F(\, ty 



(Xy)x 



dA 
dx 



X 



t 



Xy \ bt 



bXbt 



bx_ b*F(\,ty 



(Xy)x 



bXbt 



From the first of these equations 



A= 



b¥(X, t) 
bt 



\ y dy 



'b 
bt 



bx 



+xjj{x, t) 



where \jj is the symbol of an arbitrary function. 

b f bx . 



Now - - — x ^ a function of X, x, t : and y occurs in it only because it is contained 

bt H\> 



in X, therefore 



X y dij 
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r bx 




bX 


~b 


r bx 


_b 


bt. 


(Xy)x 


m 


bt 


A. 

(Xy)* 


• bt. 



bx( bx 



j Q^y) x t 



+ arbitrary function of x and t 



but the arbitrary function of x and t may be supposed included in xjj (x, t), therefore 



therefore 



But 



b¥(X } t) b 



M 



)X 



A =" L if i -s]^J ( s+^ i > 



dA_ b*~F(X, t) __d b 
dx ~ bXbt x dx bt 



M§+5 «» «> 



— — — 4-X — 
(fee &as * 6X 






+\ 



6\ 
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whence 



therefore 



b d \t d 

bt "~ dt X y dy 



h \%x[ bx 



dx~ b\bt X * [bx+ Xx b\)bt\ bX ]nf x + [bx +X *b\ 



i 



\$(x, t) 



b 3 F(x, o v be b\ x br bx , b t/ x 



6X& 



^ 



&J (x,): 



A /v ^ 



&. 



But 






f&V, rf 



a« dy\ a, 






+ arbitrary function of # and £ 



y* 



X Cv 



Before the last integration can be performed -- t~(— r") must be expressed as a 
function of X, x 9 t. If — --' be expressed as a function of X, a;, t then y can occur in it 



only through X, Therefore 



Therefore 



b 
bt 



' b\ 

(Xy)x 






6X ^(-^)=-^+*^') 



where <3> is the symbol of an arbitrary function, 
therefore 



dA / 5_ 

d%~~ \ bt 



' bm . &F(K,i)\ , X, , v , 8 f/ x 



>03 



Hence choosing the arbitrary function \jj (x, t) so that 

b 

-xjj(x, t) = &(x, t) 



dA. 
this value of — agrees with its known value. 

dx ° 

And 



ift(x 9 t)= bx<§(x, t) 






^yJt 



bxl •- W 



^/e 
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where P and Q are arbitrary functions introduced in consequence of the change in 
the order of integration. 
Therefore 



ht 



X 



•y/t 



Comparing now-— and — with their known values, it will follow that - — r 1 — = 0, 

1 ° ax dy ' dec 

therefore Q(x, t) is a function of t only and may be considered to be included 

in P(X, t\ 

SF(X, 



bt 



may also be considered to be included in it, it follows* that 



A=K(X 5 t) + \bi 



Xl r- h 



[This form of A may be obtained much more conveniently thus. 



dA dA 

Since \i-\-u\ x +v\y=0 and u——, v=— — it follows that 






therefore 



dy 



dx 



1/ 



A=K(X, £) + Sx 



% xx 



\%. 



"y/t 



The same way of obtaining A is applicable to Arts. 5 and 8. — August 30th ; 1884. 



* The form in which A appears does not appear to be related to y and —x in the same way, as would 

be expected. 

But denoting differentials of X, y, t by SX, By, St, it may be shown that 



A = M(X,t) — 



The two forms of A will agree if 



*(£)? 



Suppose that 



But 



therefore 



whence 



therefore 



AA A f A*Y 
Sxi r- f + i Sy\ r- h = a function of X, t. 



M=) 



t 



y = K, then t— =r— 



S el Xy d 



d S d S S d_^ S S ,. y _ 

cZR X y e£R X^- 



cfy \t* c&» " X it > 



y 4- $(X, £) 



so that the two forms agree. 
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But v x — u y -=f(k, i>)= some function of X, g, 

therefore 

cPA , <PA TT/X . 

where H is an arbitrary function, 
therefore 



do? dy^j 



K(\,t)+ 



dx\ 



W x 



\X 



"yh 



= H 



X 



6F(X, t) 

bx 



r bx 

(Xy)x 



The arbitrary functions K(X, t) and zr~ being implicitly contained in the integrals 

need not be expressed. 

4. Now suppose that there exists a vortex of invariable form which moves with 

arbitrary velocity along the axis of y. Let the equation to its surface be a function of 

x, and 2/— Y only, where Y is an arbitrary function of t only. 

Let 

X=L(# ? y—Y) 

then 

^ d~h 
y dy 



»<=f<-*) 



therefore 



Xi 



= -Y 



»y 



Also y— Y can be expressed as a function of X, x only from the equation 
\=Ij(x } y— Y). 

Therefore X^ can be expressed as a function of X, x only, not £. 

r bx A 

Therefore G(X, as, £)= — I — ^ does not contain t 9 since (X y )* does not contain t 9 and 

may now be written (Xy)£. 
Therefore 



(Ay)j 



SG 



S* 







In this case the equation takes the form 



therefore 



K-*+^' 



/d? . d* 



+K^)=4.- 



aH~ 



(*»£ 



6\ 



(s+*0 K < k -'>= H *• 



& 6F(\, N 
+ 



<h& 



bX 
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and now 

v=Y -A^^L. -£(k(\, t)-xt 
therefore 

is the current function. 

Of course this could have been directly deduced from the form of the current func- 

tion in the preceding article by putting r-=—Y. 

Ay 

5. To obtain similar expressions in polar coordinates. 
Let R and © be the radial and tangential velocities. 
Therefore 

u=H cos 0—® sin 0, v=R sin 6-\-% cos 

sin cos 
\ x z=z cos 0\ r \ e , X v = sin 0\ r -\- \ e 

X^ "f" ^Xtf-f- Vky = X^-f" RXy-f" _ X^ 

In this case the equation corresponding to the equation mg of Art. 3 is \ r g e —\ e g r =r. 
The auxiliary system of equations is 

d0 dr dg 

X r —\ e r 

Let differentials of the variables X, r, t when regarded as independent be denoted by 
h\, hr, bt respectively. 
Therefore 



9< 






Then since 



it follows that 



x^+RXr+T^^^' anc * sH-RflH — :ffe=o 
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therefore 






But 



therefore 



d 



™~~ T \e '{bt^btbX/ 



d 



v x —Uy=j- t (R> sin 0+® cos 0)—- ? (R cos 0-— ® sin 0) 

=?V-(Tl sin 6+% cos 0)+0,— (B sin 0+© cos 6) 



dr 

d 






— r^-r-(R cos 0— © sin 0) — #^(R cos 0— © sin 6) 



Vx—lly'- 



y dr 
dr r d6 r 



d , _v dR 



w 



dr 



dd 



To determine the current function A there are the equations 



IdA 1. /&*F(M) J>f»*r\ 



J; 



dr~~ \ e r \ b\bt bt](\ e y. 



From the first of these 



A= 



bt ~bt. 



b\ 



rbr 



+xjj{T P t) 



where xfj is an arbitrary function. 
Therefore 



rbr « d 



dK b*F(\ t) x d Br f ,„, w 



But 



whence 







<2 

dd~ 




d _ 

dt~ 


-^ + ^\ 


b 


<# X* ^ 


bt~ 


^ Xfltft? 
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therefore 



But 



rfA_ ffF(M) 

dr bXbt 

~ bXbt ' 



\ r — 



i» 



brbt 



^ v f r6r . b® f.f <rbr , o , , x 
bX -x — X r r— bX I+TW 

Ja.v r ix^J Ja.v &• v ; 



o£j 



(X#)r 



(x fl ); 



bt 



rbr , h ,. . 



OW 



or 



6^ 



*oX 



Wj 



roX~~' 



b A(X e ) 



bX\ \ $d 



rh\ 



X 



d X( d\l c^sx ^ / x^ 

dt \od0/\a J Xfl rf0\ X d 



»8 <-W/ A0 



"0 



—rr- +<£(?% <) 



X, 



where <3> is an arbitrary function. 
Therefore 

dA (&F(\ t) 



dr 



bXbt 



6 

'bt. 



l^ +r h_Q( r9 t)+~J{r, t) 



br 



Choosing Ihe arbitrary fonotion «r, so that £ «r, t) 

agrees with its known value. 
And since 



<£>(?*, <), the value of 



*ZA 

dr 



t//(r, £) = br<J>(r ? £) 



therefore 



br 



bt, 



fit 



b\ 



t 






rbr 



AA X f * AA X 






0? 



§ x p^+P(\, «)+Q(r, *) 



Wr 



A _^M + | r&r /^ + p( X? «)+Q(r, t) 



and reasoning as in Art. 3 it follows that 



A 



r 

(X, *)+] 



Xq t 



And since %.— w y is a function of X, g, it follows that 



therefore 



VV -i. jL - M/ -i. J&. * ■ "" " y r^v *"*t 



&_ , ld_ 1<P_ 

dr® r dr r* d0*. 



K(X, «)+[ 



ror I — )r 



= H 



? 1 *>\ J(**)J 



6. Now take the case of a vortex of invariable form rotating about the origin. 

Let its equation be 

X=L(r, 6—0)) 
where cj is a function of t only. 
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Then as before 



and 



b 
V 



hi' 



therefore 



CO 







^ X SK(\t) 1 d [ lz , , r 



SX 



r dO 



®=rcb — \ r 



SK(X, Q 
SX 



1 



A1 



2 



K(X, t)--d> 



and the equation in X may be expressed in the form 

d? 1 d 1 d? \l r 2 \ / 






(Mr 



6X 



The current function is 



,2 



K(X, «)~2« 



To obtain this directly from the foregoing article, put for ^ its value — c5 

7. To obtain similar expressions when the motion of every element of the fluid is 
in planes passing through the axis of z, the motion being the same in all such planes. 
Let r be the velocity away from the axis of z } and w the velocity parallel axis of z. 
The equations of motion are 

dr dr dr d I Cdp _ \ 

dt dr dz ~"~~ dr\] p J 



dw 
dt 



, dw , dw d / [dp . ~ T \ 



1 / d d d\ r dr dw 

— . | — -J-x — -4-t(J — JP I """T* — i ' == ' 

p\dt dr dzj r r dr dz 



Differentiating the first equation with regard to z, the second with regard to r and 
subtracting, it can be shown that 



therefore 



therefore 

lid 
r\dt 



'dr dvA/dw 
dr dz A dr 



d , d , d\/dw 
dt dr dzj\dr 



d , d\/dw dr 
dr dz \dr dz 



dr\ , Id , d , d\(dw dr\ 
c£s/ \^ dr ' dz \dr dz 



dz 



I d d d\ t] /dw 



n/d 
p 



dw dr 
dz 



r /dw dr\ 1 (dw dr* 
r^\dr dz) rp\dr dz 
3 c 2 



d d d 



dt dr 



dz, 
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Observing that 



This may be written 



d d d\l 
dt +T Tr +W 



dz 



r 



iy>A 



d . d . d\flfdvj dr 

dt s dr dz)\r\dr dz t 

1/dw dr 



d d d\ 

— r*T — "T" iv — )p 

dt dr dz/ r 



therefore 



r\dr dz 



d 



d 



d 



\ 



dw dr 



_L T -~--\-w-r-)< dr dz>=0 
dt dr dzj 



This becomes, on putting in Clebsch's forms, t=-~+^ ivz=^+\-~^ 



d , d , d \ \ 1 IdX dylr d-dr d\ s 

— i ^ — 4-%; — ) -< — — — 

dt dr dz [rp\dr dz dr dz 



= 



And as in Art. 2 this result may be deduced from the equation of continuity and 

, d\ . dX dX d^ dty d^r 

the two equations "^i+t— +w-r- = 0> a ~j+ * T ~1Z* W 1Z = ® on \y« 



dt dr dz 
8. Hence, supposing p constant 



dt 



dr 



dz 



r 



=■*»• *>=«*: 



+) 



^^ 



Xgz—hgr=r 



Let differentials of the variables X, r, £, when regarded as independent be denoted 
by b\, br, bt respectively, then 






{rhr b F(\ t) 



(**): 



6\ 



= G(X,r,*) + ^-° 



and 



- 1 * 

r 



i»« ~ l ~ &&X 



_ x, l pG Wx, Q 



To find the current function A there are the equations 



1 /fe 2 F(x, Q 



6 f r£>r \ 

bt) (\\5) 



\_1 / 6 2 F(X, Q 
X 2 r \ &X6i5 






1 dA 

r dz 

1 dA 
r dr 
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From the first 



therefore 



But 



whence 



therefore 



But 



therefore 






bt bt, 



frbr 



J 0w) 



dA 6 3 F(X, t) , d b 



dr b\bt 



X r ' 



dr bt 



f 



b\ 



rbr , d , , . 



(**) 



$r 



* 






6 
ir 



+ ^ix 



A. 



^ 



+X 



*&x 

b 



t 



b\ 



d \*d 



bt dt \? dz 



dA_ b*F(X, t) _ bCrbX 
dr~ bXbt r bt](x z ) 



bXbt Xr btlfrt Xr bt, 



'rbr b ,, x 



b_ 

bt, 



' rbX 
(X z )r 



4 



rb\ 



bf 1 



W 



V t 






rbX d I X t 



X» dz^ 



dA fb*F(\t) bfrbr\ X t b 



Hence, choosing t|»(r, tf) so that 



< Km)=t-# , >0 



£>r 



and therefore 



i//(r, £)== k>(r, £)br= rbr 



As/* 



br 



bt 



rbX 



Ow) 



y ,i 



rSr( ^);+IW 



Pw)r 



+P(X, *)+Q(r, 



flJA . 



the required form for — is obtained. 



dr 



A: 






+ W^);+P(^ 0+Q(n 
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and reasoning as in Art. 3 it follows that 

A=K(X, t) + 

ut - ( 7- is a function 01 a, g. 

r\dr dzj 



ror( — y 



mi f* 

I nereiore 



therefore 



r\dr % r dr dz 2 



rAdr % r dr ! d$J 1 x ; J \ W* 



H[X, £ 



H 



X,- 



b¥(\ t) 
o\ 



rbr " 



9. For a vortex of invariable form which moves parallel to the axis of z. 






where Z is a function of t only. 
As before 



^=-Z and ^=0 



T= ^!%i) = iA( K(M ) 



4^. 



,2 . 



r 



8\ 



r d 



z 



r 8X r dr x x f 



2 



2 . 



Z 



Therefore the equation in X becomes 



1 d 



,2 , 



bt (/• j 



;0 K ^^-^ Z 



'1 d_\ 8 1 1 — 1 



+^ (k(x,*)-*2-z =h(x, -J^ + 



rSr 6F(\, *)' 



(Xs), 



b\ 



n2 » 



K(X, «)-s-Z =H \, 



6X 






therefore 



i/# . * _i d\ (K( ))==H / x **&,«) 



& 



rof 



(K) 



The current function is 






To obtain it directly from the preceding article, it would only have been necessary 
to put 

Z 



x* 
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Illustrations. 
10. Example I. Take the simplest case of the equation in X given in Art. 4, viz., 



where c is a constant. 



\dx 2 <% 3 / 



/oft (y—Yy\ 

It is required to examine whether \==(y , )./—+ — ^~ / can represent vortex sheets 

in a motion, part of which is rotational and part irrotational (f, a, b being constants). 
Substituting in the equation 



Also 



J a 1 



To find \ji it is necessary to solve the equation 



*tt _l f 2 (y~ Y ) ^i(y— f 2 ~\^=Q 



The auxiliary system of equations is 



One integral of which is 



and the other 



T ~ 2/(2/ - Y) ~^~2A~ T 
5 3 a 2 



(/)(~!+^? : )= const. =m 



ab . , a? f'f 

t— 777 sm x - a/ £-= const. =w 



where for m must be substituted its value. 
Hence 



^(f^+^ : ) 



oc 



a& . , a 

\fj=:t— — ~sm * 



V: 



2/ , A 8 + (y-Y) s 



a 3 & 3 
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Whence 

a?~T~ b* ') dx~~' 2f 
+ 



a? ' ¥ 



J dy 2/ 



Substituting for u and v in the dynamical equations 
■-J-V=-^ (^ 2 +(?y"~"Y) 2 )"~^Y4- an arbitrary function of £ which need not be considered, 

Now the equation determining x which is 



dt 



becomes 



^ x 2/(y-Y) ^ x fa_2Mdx 

The integrals of the auxiliary system are 

(/)^$+ V: V/ ==conBt - =m 

ah . _-.# // 
^ — 7-isin x - A/ ~= const. =w 



2/ a y m 



( n\ f( —-"J) 4# f r nv * 2/ 

X== T^J" sin ^('" w ) + 6 V / Y cos ^(*-^)+ a Unction of *. 

Hence one value of the integral of the equation in ^, which may be called x> is 
obtained by substituting for m and n their values, and is 

=^>(y--Y)+*(y-Y) 

W hence 
But 
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If now e be some constant, then 



u 



5+68 ! J^y-X y= -^^/-^"J + ^ 



2/' 2 (a 2 + & 2 ) 
Choosing the constant e= ~rz , the right hand sides of these equations 

become e\jt Xi e\\s y respectively. Therefore 

dy\ X J ' a?W V J aW \a*~ r b* Jdy 



Putting in for x' an( l $ their values, it will be seen that for the x 0I * Clebsch's 
forms of expression for the components of the velocity, it is necessary to take 



x 



a y-/%-Y)+Y(2/-Y)-2/^^-|sin^ - 7^~fy 2 



v 7 



a 2 Z> 2 



The terms, containing £ only, may be omitted. 

Thus if p be the density of a distribution of matter of which x is the potential 

a 4 — ¥) x(y—Y) 



H 4nr\d& T dfj 4tt 



aW /x^ (y-Y) 8 \ 8 

therefore 

c a 3 — Z> 2 «(«/ — Y) 

The potential of this density (see Art. 15) at an internal point is 



L-, 



X 



c a—b , TTX ca& i . _, a , a? 

■ #(y — Y) + — sin x /-s— 7 — -^ 2 — sin"" 1 



2a + & M ' v * 1/r 2r A 2 (y- Y ) 2 ° 1XX y/^-j-^-Y) 2 

\ V a 2 + Z> 2 
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Adding to this the cyclic term 



(jIXiO » -i Jb 

— sui l 



2 y^ + fy-Y)* 

giving a cyclic constant --7ra&c=7ra&(+2£) where 2£=~ — —J the expression 

obtained is 

x 

T^i "to- Y >+ T 8m 773T5=W 



V 



a 3 s 3 



This will not agree with the value of x ? unless 

c a— 5 & 3 — 5 3 . 

2 # a + & — a 3 5 3 ^ 

ca&_ /(a 3 + Z> 3 ) 
2 ab 

Y=0 

The first and second equation require a =6. 

Hence this method will not lead to a determination of the irrotational motion out- 
side the cylinder. It does not prove that there is no irrotational motion outside 
continuous with the rotational motion inside the cylinder. 

Supposiug Helmholtz's method applied to this case, it would be necessary to find 

G 

a value of A which is the potential of a distribution of matter of density — — inside 
the surface 



^ . (y-Y) a 1 



The result is that inside, 

and outside, 
cab 



cab 



a(a + b) b(a + b)^ 



A ' 2 



c+ HW^ + vw^ )+ ^^-^§^(£ r ^ 



where 



#2 C^ — Y) 3 
-r — -+~rz — — = 1, and C and C are constants. 



If the constants C and C be properly determined, these expressions will be 

$ 3 (y-~-YY 
continuous at the surface ~+^ Tr x =l. Their differential coefficients are also 

continuous. 



<z 3 
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_ T dA ca(y—Y) dA. —cbx . . .. ., ,. , 

JNow — -=■— — — - — = inside the cylinder. 

dy a + o dx a-\-o J 

dA dA 

At the surface of the cylinder. — and -7- have the same values whether calculated 

d dx ay 

from the value of A inside it or outside it. 

In order to have the rotational motion continuous with the irrotational motion, it is 
necessary that all over the surface of the cylinder 

•*' W a + h 

• 2fx cbx 



a? a + b' 



But these equations cannot be satisfied unless a=b, Y=0. 

The solution may, however, be completed for a finite portion of the plane of x, y 
outside the cylinder by means of Example III. which follows :— 

In Example III., put c&=0, this will make 

X=A of Example IIT.=/.~~ log (v^+c+ s/b + e)+J^,(x*~(y-~Yy) 



where 






x * _ [ >-Y) 2 _ 1 



2^ 



a 2 + e W 1 + e 



Therefore the current function A of this example is 

k-xt=J^ log (^¥+l+v , ¥+J)+^L,(o?-(y-Y)<>) 

This is equivalent to the form in the abstract printed in the Proceedings. 

For the motion to be possible, it must be supposed to be confined to a cylinder of 
finite section, appropriate surface conditions being supplied at the surface of the 
bounding cylinder. 

11. Example II. In the preceding example it was shown that none of the methods 
would apply for the whole of space surrounding the rotationally moving liquid. 
Knowing Kirchhoff's investigation of the rotating elliptic vortex cylinder, let the 
form of the equation given in Art. 6 be considered, and the following simple case of 
it taken. 

3 D 2 
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A particular integral is 

_ , ^ J> 3 cos 3 (6 — co) r 9, sin 3 (9 — co)T 
A,_(/)| - + — J 

if 2/(~^+v 2 ) — 2c5=c(= — 2^ of Kirchhoff's " Vorlesungen liber Mathematische 

Physik." Zwanzigste Vorlesung). Therefore cb is constant. 

Let x, y' be the coordinates of the point x, y if referred to the principal axes of the 
ellipse which rotates with uniform angular velocity cb ; v! v' the components of the 
velocity parallel to these moving axes, 

Thus 

x=x cos cbt-\-y sin cot ; y=- —x sin cot+y cos cot 

The velocities along and perpendicular to the radius vector are 

R =; !( x -t)= - </>(H> *> *»— > 

^=Rcos (0— o)) — ©sin(#— <y) = ( — -— c5 W 

*/=R sin (#-»o/)+© cos (^— «)= — (J— fl>V 

w=w' cos cbt—v sin a>£ 
-y=w' sin d>£+?/ cos cot 

a'* , ?/' 3> 



*=(/>br+ 



2 



a 2 b 



To find xjj it is necessary to solve the equation 

cl^r dyfr ^ _ A 
at ax ay 

The auxiliary system is 

dt dx dy d^r 
1 u v 

It will be more convenient to obtain xfj in terms of x\ y\ 
Since 

dx — cos cbt.dx-\- sin cbt.dy-\-cby' \dt 

dy= — sip. a>£.c£x+ cos cbt.dy—cbx\dt n 
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therefore 



therefore 



therefore 



di 


dx f 


1~ 


Vj cos cot +v sin cot -f <%' 


dt 


dx dy' d\jr 


1 


u + ooy f v' — cox ; 


dt 


dx f dy f dijr 


1" 


y'~~~ a?'"*"" 

4/^3 ""4/^2 



dy dtyr 



•^sin d>t + v cos c5£ — d>a/ 



The integrals are 



therefore 



C/)(5+«f)=» 



, a& , (x f / f\ 

2/ \a V m/ 



X =(/)(^+-yT 



7& 



, L ab . 
\jj=t — ^.sm 



-l 



xf 



a 






The current function A=X— ^cbr* 



= (/)(S+S-f(^ + ^) 



# 



To find ^~+"V" 

(d d d\ 
dt^ u d~^ V d^ ) * n a ^ orm * n w ^^ x > V'> t are ^dependent variables, 

A +u A +v l\=Z. + ^ AjlM. A +u ( d A A + W A\ +V (M Ajl^I A 

K dt dx dyj dt dt dx' dt dy'* \dx dx f * dx dy'j~ \dy dx~ dy dy' ) 

7 7 

= j +(coy'+ic cos cbt+ vsiri d>t) ,^(~<bocf-~u sin cbt-\-v cos cot) — 
dt K J ; dx l K ' 'dy' 

—A i W. A-J^A A 

dt & 3 dx a dy 

therefore the equations 

du du die d (jp ^ 

at dx dy dx\p 

dv . dv , dv d fp , __ 

<fe dx dy dy\p 
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become 

d 2fy' d 2fw' d \ , , . x d fp 



di W dx a 2 dy 



~)(u' cos cbt—v' sin <££)= ~~ 7p+V 



d , 2fy f d 2 fa/ d \ , , . . , , . x d (p , _ r \ 

-+~T7""T7 — tttP smcot + v COSO)t)=— — -+v) 
v &£ o 3 cfar a z dy ] x ! dy\p J 

therefore 

• ., ,/ 4 / 9 4 /^ i -a\ i •* '/ 4 / 3 . 4 ^ -a\ d /P lV\ 



cos ^y^_ 4 ^ +l fi2)+ sin &t.af (_g+^-c6»)= -^+V 



therefore 



\ ercr er / <x<2?\p / <%\P / oas \p 

therefore 

■?4-Yz=x /s (|c5 2 — 2 +^)+2/ /2 ( i^ 2 ~""i^ — h ix3/+ an arbitrary function of the time. 



To find x 



p ' 2V ' cfe cte (% 



The auxiliary system of equations to find x * s 



dt dx dy d^ 






jp 



Substituting for u, v, — +V their values, there may be substituted for these the 
equations 



One integral is 



/— - J- -f" - am 



Another is 

2? 



\a V m J 
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To find the third, substituting for y' its value b A/^— ^j 

J 

dx' ___ d% 

therefore 

therefore one value of x satisfying the partial differential equation, which may be 
called x ' m 

/y2 7)2 

/ /» M[ ^_ / / 



Also 



/ T '2 n/% 

. ah . , <z 



Whence 



2/"*~ A 72 .'V 



a? W 1 j dx f 2f 

V 2 y' 2 \dy]r a/ 



Thus 






If, therefore, 6 be put =2/(c5— ■■— - 






Thus the proper value to take for the x of Clebsch's forms is x'+#- Omitting, 
as unnecessary, terms containing t only 

cft — h'* I f f\ ~a 



a3&3 -^ — y~ a 2 jay— /^2 y /3 



V 



a 3 5 2 
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Calculating 

p 4tt\ dx* "*" rfy V 47r\dte /3 " r ^' 3 , 



c / 1 1 \ a/y' 



47r\& 3 a*](^ i£\* 

a? b % 

The value of the potential due to this density is the single-valued expression given 

in Art. 15, viz. : — 

(x f 
. ^ a j a/ 

sm 7^77 ~ sm ~ 7^+F 

Add to this the term — sin~ ] --/====| to give the cyclic constant — (— 27r)=7ra&.2£ 

in Kirchhoff's notation ; then, in order that what is now obtained may be the same 

as x> it is necessary that 

/.& 3 ~& 3 __£ a—-b 
* aW~~2a + b 

the latter of which is known to be true, but the former will not also be true unless 

•__- 3/_ cab _ y ®b 
~~a£~~~~(a + by = ^(a + by 

Up to this point no relation has been assumed between c5 and f. 

Supposing, however, this relation satisfied, the velocity potential at an external 

coin hp 

point is obtained by adding the same term — sin" 1 ,-% — ~| to the potential found in 

Art. 15 for an external point. 

Thus velocity potential at an external point is 

cabx'y f I & 2 -f£ 3 + 2e \ , cab . , x f 



(o=—~=z — , ,^=2^^ , 7 . 2 in Kirchhoff's notation. 



^-& 3 \2 v /(a 3 + e)(& 3 + e) P 2 yV + e 



x' 2 y' 2 

a? + e W + e 



where e is a root of the equation ^37;+^^== 1? the axes #', ;/ turning round with 
angular velocity g>= — 



(a + bf 

The expressions for the velocity, which may be deduced from this, might also have 
been obtained by Helmholtz's Method from the current function A which is the 



OF WHICH IS MOVING ROTATIONALLY AND PABT IRKOTATIONALLT. 393 

potential of the density — -— throughout the cylinder 0+^ =1, and which is given 

in Example I. 

v 
These values make-+V continuous at the surface. 

P 
12. Example III. Supposing the relation found in the last article between/* and d> 

not satisfied, it is required to find a solution if possible. 

Using the elliptic coordinates e, v which satisfy the equations 

+ 7f~~ =: l 



where 



Hence 



and putting 



a 2 4 6 b 2 + e 
a 2 + v b 2 + v 

— 6 2 <€<00 

— a 3 <v< — ¥ 

(a 2 + e)(a 2 + v) (b 2 + e)(-b 2 -v) 






X ~ ~ ^-6« ~~ ^ ~ a 2 -b 2 



a= log (\/a 2 +e+ v/& 2 +€) 
/3= — tan" 1 a/ 



— 6 3 — t> 
+ <x 2 + v 

cPY d 9 Y 
Laplace's equation becomes -^ + ~— =0. 

1 da z d/3* 

If V be a function of e only, then V=Ca+C Similarly V=Cj8+C / is a solution, 
But if V=E.U where E is a function of e only, U a function of v only, then 

d^ JJ3P 

Suppose there exist a value of E such that -r^==p s E. 

Then 

d 2 U 

Therefore if 

E=A^ a +Be""^ a 

U= A' cosjp/3+B' sin jp/3 

E=A(v / a 2 +€+ v / 6 2 +€)^+ B(n/^+c+ v/S^+c)"^ 
U= A' cos s^> tan" 1 a/ ^-^ > — -B' sin -4 jp tan" 1 a/ 
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When 

p—1, E=C\/a 2 +e+D'v/6 3 +e 



V=C\/d 2 +v+T>'V-P-v 

U = G'( v+ c ^~\ + D' v/"(a s +v)(-F^) 

and so on. 

Now suppose it required to express A as the current function of an irrotational 
motion under the condition that 



^=f-,+r^-2( x " 2 +y'") ^erever -+--: 



Assume if possible 



f x i+f y ¥~i v+w 



A+Bhg(s/a i +e+V / b z +e)+ Clj+^j+D^+eK^t) 



, <x 3 + 6 3N 



The form of the expression on the right-hand side shows that it is the current 
function of an irrotational motion. 
But 

cfi + b 2 ^aP — b^ / x ,% y^ 

therefore the assumption is 

= A+B log (v / ^H+ v'&^R +C ^ (a?*-^) -C fc 5 ^ 

+D — - — \/(a 3 + € )(^ 2 + e ) 



k a 3 -f-e & 2 -fe. 

If this is to be satisfied when —+^~=1, and therefore when e=0, it is possible to 
add Jc f — +J7~- 1 ) to either side and then equate the coefficients of a/ 2 , 2/' 2 . The con- 
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stant A can always be chosen so as to make the absolute terms equal. . The elimination 
of h from the two equations for 0, D in which it occurs will leave the one relation 
between them which must be satisfied. 

/V 3 ip \ 
Adding & (~y+T£""-"l ) to the right-hand side, the conditions are 



f d> ^at-V* , ^ 7 a 2 ~& 3 , k 

- — .-==C — — h-D.a&.-— — ^-~ 

a* 2, 2 2a 3 a? 

j2 2 — °' 2 2& 3 ~" + "& 2 

and an equation to determine A. 

The velocity parallel axis of x r inside cylinder =-77 = (--T—<b J y'. 

The velocity parallel axis of y' outside cylinder = — — = — •( ~^— co) a; . 
From the assumed value of A for space outside cylinder 



^ = _C(a*-&%'-D(a*-&%' ,sj 



dy >— ~\~ -/j -v- -/* ^ V F+^ 



1 "* fr.i.-n^-^V ^ 4- y " 



& 3 + e 






^ fE-l-D^^Vf + y " 



For the continuity of the values of —„ l —~, at the surface e=0, it is evident that the 
coefficients of — , — must vanish. Therefore 



B + D^ = 



Also 



|f— c5= -C(a 3 -6 3 )-D | (« 3 -& a ) 

^f— a>= C(a a -6 3 )+D-(a 3 -6 2 ) 

3 E 2 
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These are the same equations as those obtained for the continuity of A if &=0, 
whence 



CvvLf 



(a 2 — Iff 



4/_ . a? + W 
ah ah 



ah f la? 1 -4~ o^ 
hererore 



r '2 ? /2 A 



const. +a&|/~ a |f-c5j log (\/a a +eHVP+e)+i.-^^ -c5^~--)(a' 3 -?/' 3 ) 



a& / ,71 , 1\ \ /7 3, , wis , \ Y ^ y' 3 



/ ^+i5 -^ Jv'V+eX&He) 



a 2 — 6 2 \ \a? W 1 } l ' '\a 2 + e' & 3 + e, 



To examine whether this will make the value of --+V continuous at the surface, 
Inside it is known that 

<n / 9 fA 9 /2\ / 9 /A 9 /2\ 

^+V=a;' 3 (^o> 3 — :- r +_J+y'a (i« 3 &^+^¥/ + arbitrai 7 Unction of t. 



Outside 



*^=-f -*((£) '+(*)> «*»* *»*» ° f «■ 



where <£ is the velocity potential. 

Now as a> is constant, if <jf> be expressed as a function of x 7 , y (e is a function of x\ y) 
then t can only occur in <f> through occurring in x\ y, 
Therefore 



d(j> __dcj) d%' d$> dy f ^ , . - , , . , 

-_ ~ _— — - * --—- —J— ~~~~~. * —— w *0)iJ ""-*" V • 0) 00 

at dm at * ay dt u 



<y 

But it has already been shown that the velocities are continuous at the surface* 
Therefore at the surface 

^ — >(L i -\ d $ — o '/Z i - % 

Therefore at the surface 

=+ V= — 2a/ 3 / ^|— Jc5 s J — 2y'% (4j— i^> 2 )+ arbitrary function of i. 
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The difference of the values of -+V" as given by the above expressions obtained 
from the motion inside and outside is 

/ 1 1 \ 1 /x' 2 ?/ 3 \ 

2 ^ V + &7 ^J \a* + W) + arbitrar y function of «• 

As— +78 = 1, choosing the difference of the arbitrary functions of t to be constant 

The current function of the irrotational motion is A, 

Therefore the equation to a surface always containing the same particles is 

# 

X = - (x^ + y >% ) + A = const. 
That is, in elliptic coordinates : — 

-(c+v+a a +& a )+a& \f — —-co I log (v / ot 2 + c +v /r & 2 +c) 
ah 






(a?-wy 



-4 — co — —€+—-—--. 2( f —^r— ™ )v (a +€)(6- 3 +c) v+ — — )= const. 



.a& db \ x 2 V a?W 






Now if e be chosen so as to make the coefficient of v vanish, it will always be 
possible to choose the constant so that the equation is satisfied. Therefore the 
elliptic cylinders corresponding to the values of e which make the coefficient of 
v vanish will be parts of surfaces X=const., and will therefore always contain the 
same particles. 

These values of e satisfy the equation 

Solving this in the ordinary way, the roots obtained are 

€=0, € = 






It is necessary to examine whether these roots satisfy the above equation, or the 
equation obtained by changing the sign of the radical on the right-hand side. 



398 PEOFESSOE M. J. M. HILL OJST THE MOTION" OF FLUID, PAET 

As positive values of e only are required, it is necessary to examine the relations 
between J> and f which will make e positive 

\ a% % f 



co 2\/co 2 
+ 



/ * ccbj\f ab, 



If 



d) 2 .... 

oo < ~ < — ~ € is positive. 



f a ^ 



e l 



is infinite. 



2 co 2 . . . 

— <— < -r € is negative. 

ab j ab ° 



e is infinite. 





• 

CO 


2 
ab 


2 
ab 


CO 

f 


a? + b* 
aW 




CO 


a? + W 



e is positive. 



/ & 3 & 3 



e is zero. 



& 3 + & 3 w . . . 

— — - < - < oo e is negative. 

a 3 6 3 / - 



If the relation between co and / is such as to make e positive, it is necessary to see 
whether the equation in € above is satisfied. 

ia? + & 3 co 



.// a3 fr 3 



(5 2\» 
a 1 

,2_L \ / «/ 



a 3 4"€=-T-r 



co 



> 3 4 



// a?b 2 

/,© 2\» 

& 
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First suppose 



Then 



6) h z 

CO <-;<—--- 

j ab 



6) 2 

— a — {- 



V or + e= /t^TT^ a 



f) a%^ 



-&*+* 



V 7 "^ 



so that 



y/ a a 6 3 



, — a— + -|( — o — t-- 



\fj a?W 

and the equation in e is satisfied 
Next suppose 



Now 



2 d> a? + 6 3 



ca 2 

a 






b* + 2 - 






\fl a ^ 

These values, however, do not satisfy the equation in e. 
Hence it appears that if 

& 2 

—.00 <-7< r 

/ cib 

the elliptic cylinder for which 

\ a?V* f 



v// ^ 3 
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is part of one of the surfaces X= const., and therefore always contains the same 
particles. 



Hence, if the smooth hollow rigid cylinder 
uniform angular velocity c5, such that 



x 



a* 



■ —mi— 



y 



&*+. 



1 rotate about its axis with. 



41 



€: 






7 ' 



4 



-, and 



d) 2 

oo <! — <** — - — . 

V^ aV 



a% % 



x f% . y f< * 



it is possible that the fluid inside the geometrical surface ^f +t^ = 1 should move 

rotationally, and that the fluid between the two cylinders should be moving irrota- 
tionally, the rotational and irrotational motion being continuous. 

The components of the velocity of the rotational motion parallel to the axes of the 
sections of the cylinder by the plane of x, y are 

,/2/ .\ , ,/2/ 



The components of the velocity of the irrotational motion in the same directions are 

dA dA 
dif dx f 



t'.t'., 



aby f 



and 



a 3 --& 3 



abx' 



2 {^ + ^)- <b } VwVe~ 



a? — ffi 



*{f\j+v)- & } VS~ 



4/ . a 3 -f & 3 

ab ab 



4£ .a 3 + Z> 3 ' 



These expressions for the velocity will not agree with those obtained by Helmholtz's 
Method, viz., u^^rJ \/t^-l) and v=^J yy/^-l), unless «=-^, 



W 2 \'V 6 3 -fe "7 a 3 -& 3 \'V ^ 3 + € 

i.e., in the case of Example IT., and the irrotational motion as obtained by Helmholtz's 
Method is not continuous with the rotational motion. The other method suggested 

1 Id / v d "v\ 

in this paper (in which the potential of the density — — ( -~ +vf ) is calculated) does 

not lead to a result. 

13. Example IV. To consider the case in which the vortex sheets are coaxial 
circular cylinders, and the molecular rotation is a function of the distance from the 
axis. 

This investigation will illustrate the reduction of the components of the velocity to 
Clebsch's forms. 
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In this case ~r^-\ — -; — I — g ~Jm := a.function of r only, and X=F(r). 

It is supposed here that the initial line from which 6 is measured is fixed. 
Let R, © be the radial and tangential velocities, so that 

r dv 



®= -!=-*» 



xp satisfies the equation 



t+ R f+®s=° 



The auxiliary system of equations is 



dt dr dd dty 



One integral is r=m, the other is t-\~ , 0=n 
Therefore the integrals are \=~F(r), \jj=t+=^-~0 

J: \T) 

The most general form of the integral will be ^=^(^1, n) where for m and n their 
values must be substituted. 

Substituting for R and ® in the dynamical equations 



V 
P 



+Y= (F(r)} 



2 



dr 



r 



To find x> 



% + B £+• S=« w-!-^ 



Therefore 



f+°l-™2=^w-f(F«) 



2^: 



The auxiliary system is therefore 



dt dr dd d% 
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«FW} 8 -|{F(r)} 8 7 
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The integrals are 



r=m, 



r 



Jb (r) 



X 



-0 



'-^^^'H- 



therefore the integral of the original system is 



X= 



( , ~ ^ — ftP-J - *— — — " 



F'(t-)J 



{F(r)} a -[+X(m,w) 



Take X=F(r), xfj—^m,, n) 

To reduce the components of the velocity to Clebsch's forms, it is necessary to find 

X and "9 so that 

rfr f dr 9 rd6 rcffl 



t/*G* 



Whence 



0=e ^_^ +{ n^ )lmw 



dr 

r 



JL 






bX nr)-rWXr ) () l 
^bn {F'(r)}» T W ] 






6m 6m 



"F'M 1 



F(r) 



( FM 1 »- -I- — — 4- b - - F ^ 

1 v/i r ^6w'F(r)^ bn¥(r) 



bX , _. .6^ rt 



&m 



6m 



bX 



+F(r)~=-i(FW} 2 -f{F'W} 2 7 



To satisfy these put W=n®(r) in each, then since m=r, it follows from the first of 
these equations that 

X=-n\¥(r)&(r)dr+®(n) 



Substitute in the second equation, therefore, 



F(r)0>'(r)dr+e'(n)+F(r)*(r)= -i{F(r)} 3 - {F(r)} 3 T 



It is necessary that %\n) should be constant, and it will be found that the constant 
may be taken &s zero. 
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Differentiate the last equation with regard to r. 
Therefore 






r 



Therefore 



X 



=e 



nr)-rnr)V™ 



+* F(r)F'(r)4{f'(r)f+ 



F(r)W(r) 



r 



{W(r)/j 



X=F(r) 



^=-(^+F^)( F »+ ] T : 



If the rotationally-moving liquid be bounded by the cylinder A.= 0, and its radius 
be r=a ; then F (a) = 
Therefore when 



r=a, x =0(-aF'(a))+t\ -^'(a)) 2 - (F(a)) 



k 2 



da 



a 



A suitable value for the velocity potential at points outside the cylinder r=a is 



^=(?(-aF(a))+«j-i(F(a))»-f(F(a)) 



2 



da 



a 



This will make the velocity and pressure continuous at the surface of the rotation- 
ally-moving liquid. Also the velocity at infinity will be infinitely small. 

14. Example V. This case is of interest, because one set of the vortex sheets, viz., 
ar 2 (z-— Z) 2 +&(r 2 — -a 2 ) 2 =const., consists in part of ring-shaped surfaces. The results 
only are given. 

If a and b are positive, and the constant <6a 4 , then this represents ring-shaped 
surfaces. 

The equation in X of Art. 9, includes as a special case 



lid? Id, d\, 
r^yir* r ar dz /, i 



A particular integral is 



\=ar 2 (z~- Z) 2 +6(r 2 — a 2 ) 



2 



giving 



2ar(z 



Z) and ^^Z~2a(^-Z) 2 -46(r 2 -a a ) 
3 F 2 
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The value of ~+V is 

P 

where 0\t) is an arbitrary function of t. 

The differential equation in xfi has two independent integrals one of which is 

ar~(z— Z) 2 +&(r 2 — a 2 ) 2 — const.—e 
and the other 



V«\/«-^-* 8 ) 8 "" const -— / 



where after the integration is performed, the value of e must be substituted for it. 
The differential equation in ^ when solved will give 

x ^ ar 2(^Z)~fa(^Z) s +4&a 2 ^ 



46 s + ah 



r*dr 



-2&(4Ha)«H- 7r J 7nF ^ l +G(,/) 

where G and <f> are the symbols of arbitrary functions. 

Finally, in order to express r as -j^+\-~^ and w as ^^\-j~, G(e 9 f) is taken as 

(8b+2a)(b^--e )f. 

Then X= (8& + 2a)(e — e), */>—/, x == a ^ ove expression with the value of G(e,f) taken 
as (8b-{-2a)(boL 4f ~-e )f; where 6=ar 2 (^—Z) 2 +6(r 2 — a 2 ) 2 and 



, 2^/a\/e— b(r 2 — a?f* 



^ is therefore 

a^( 2 -Z)-|a(,~Z)H46 a ^-Z)+^+ 4&! 7rf77=S 



<y/# J v^ 6 -~ ^ S — °^)^ 



'b* A — e f 



dt 



— (ib4-a){ 7=^ ) ,-• r T-^^-^^r an arbitrary function of t. 

The value of — — ( ~+ — , + vf ) calculated from this is complicated. The writer 

4:7r\dr z r dr dz~ J x 

has not succeeded in applying any of the methods of this paper to complete this 
example. To complete the solution it would be necessary to find a value of the 
velocity potential <f> which is continuous with ^ all over the surface e= e , and then to 
examine whether the rate of variation of <j> is equal to the rate of variation of ^ 
normal to surface e=e . The former part of the work is always theoretically possible, 
but it may happen that the latter is not. 
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The values of the components of the velocity found however completely solve the 
following problem : — 

A hollow, smooth, rigid surface of annular form, whose equation is 

ar 2 (z— Z) 2 +6(r 2 — a 2 ) 2 = const. 

moves parallel to the axis of z with arbitrary velocity Z, to find a possible rotational 
motion Lide it. 



Appendix. 



$ 3 y* 



15. The density inside the elliptic cylinder —+7^=1 is 



a? 6 3 

Gil 1 \ 

— \ X y 

__ 4.tt\W a* J J 

it is required (on account of Examples L, IL, and III.) to calculate the potential inside 
and outside. 

(It may be noticed that although the density is very great near the axis of the 

elliptic cylinder, yet the total mass of matter inside the cylinder -+-~= const., 

(aj 

however small the constant may be, vanishes. Hence it is not singular that the 
potential should be finite). 

The density varies as xy on every ellipse whose equation is l -+™ = const. 

It will be well to commence by finding the potential of a cylindric shell bounded by 
the cylinders 

where m is a little > 1, and where the density varies as xy. 
The following are suitable values for the potential 

a— 8 
V'=:C.A — ~~zxy inside 

and 

tt r* a ft (t a 3 -f/3 2 + 2e \ ,. - 

V =(3.-zr—T[ 1 — - /t ~e — rr^ — = )xy outside 
where 

^ 4- y * = 1 
cP + e^QP + e 

and C has to be suitably determined. 



4.06 PROFESSOR M. J. M. HILL ON THE MOTION OF FLUID, PART 

For V is finite and continuous inside the cylinder and satisfies Laplace's equation . 
It is continuous with V at the surface of the cylinder. 

Also V satisfies Laplace's equation outside the cylinder, and — , — both vanish 

at infinity. 

To find the volume density p of the shell, let Sn be the thickness of the shell, dn 
an element of normal drawn outwards. 

Also 

dn dn ' 

Calculating the value of — at the surface l -z +75= l,i.e. 9 at the surface ~ — 4-^o = 1 

6 dn a 3 /3 3 a 3 -}-e /3 2 + € 

where e=0 

dY n ot—{3 xy 



rf^ • r ° '* + /3'*ft 9 



x % 7/ 3 



£> being the perpendicular from the centre on the tangent to — +— == 1 at the point x, y. 
Similarly at the same point 

therefore 

dY dY' pGxy a 3 — /3 3 
dn dn 2 * a 3 /3 3 

Also — —— where a+Sa is semi-major axis of external boundary of shell. 
Hence 

Gxy a 3 — /3 3 a 



P— 



• iv • 



87T a 3 /3 3 S# 



Now consider the cylindric shell bounded by the two surfaces 



JL, 3/ 2 =land ^ , 3/ 3 =1 

(m&) 3 (m&) 2 [(m + Sm)&] 3 [(m + Sm)&] 3 



The density at the point x 9 y inside it is 



1 1 

c \& 3 a? y 
& 3 5 3 
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To find the potential of the shell, put 



where 



a? ' W 1 



a=ma 9 B=mb, Sa=aSm, ~+t^=^ 2 



Hence 





C= — 2c— - 
m 6 


Therefore potential inside is 






2cBm , a — b 
mmmm — .jl. npii 

m 3 2 a -f & ^ 


and potential outside is 




2c8m ab 


/ m 3 ^ 3 + m 






m s &2_ a 2 ^ 2 v /(m 3 a 3 + €)(m 3 & 3 + e) 



where e is given by the equation 






mhfi + e w 2 6 2 + e 



To calculate the potential of the whole cylinder, put e=m 2 P so that 



•^ , f ___2 



+irrs=w 



a 2 + P ' & 3 + P 

Then the potential of the shell, considered above, at an external point, becomes 

__2cSm t ab I __ a 3 + 5 3 + 2P \ 
m 3 V~a 3 \ 2 v /(a 3 + P)(& 2 +P)/ ^ 

The limits of integration for m are and 1, when integrating to find the potential 
of the whole cylinder at an external point. 

Therefore the limits of P are oo and that value of e which satisfies the equation 

/v>2 qi" 

„_l^ — = i an( j which makes both a?-\-e and & 3 +e positive. 



Also 



(<x 2 + P) 3 ' (& 3 + P) 2 
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Therefore the potential for the whole cylinder at an external point is 



ab f e (^ 

00 \ 



X* 






r 



a 3 + 6 3 + 2P \ |> 3 + P) 3 ' (6 3 + P) 3 



2/(fl» + PXJ» + P) ( 



«r 



" «■«■&«•* 



t 



a a + P ' & 3 + P 



2 



-cEP 



cab 
¥^a? Xy 



1- 



g 3 + 5 3 + 2P 
2 V /V + P)(5 3 + P) m& tan _ t 



05 



2 



+ 



r 



2 



< - 



V 



5 3 +P 
a 3 + P 



cab 



ffi—a* 



xyll 



» 3 -bP & 3 + P 



a 3 + & 3 + 2e 



00 



£C 



, ,-f-ftaiH 



V 



& 3 +e 
<x 3 -f-e 



tan"* 1 - 

y* 



cab I . o? + b 2 + 2e 

xy 1 — 



& 3 -~ a 3 



, cab / . -, 
+ — ( sin"" 1 



03 



2v / (^ + e)(^ + 6)/ T 2 \ -x/^ + e 



sin l 



iXj 



s/x^ -f £/ 3 , 



«2? 3 ?/ 



To find the potential at an internal point x, y. 

Suppose that this point lies on the elliptic cylinder 2 ( 

Then the potential at as, y = potential of matter inside this cylinder 

+ potential of matter outside it, 



The potential of the matter inside 



x 



2 



r 



(jiaf (fib) 



1 is obtained by taking the same 



integral as in finding the potential of the cylinder —-{-~=l at an external point, but 

Lv (J 

the limits for m are now and //,, i.e., P= oo to P~a root of the equation 



— — -./xr 



a 8 + P"*> + P 

which makes a 3 +I\ 6 3 +P both positive, but this root is zero, since 7-7^+7^=1. 



This gives 



W ' W 



cab xy I (a-~hf\ cab/ ^bx _^ 



t/.&.y 



c a—b xy , cab 

2a+6/* 2T 2 ^ 



sm 



-i 



OS 

a 



v 



x 2 3/ 3 
a 3 & 3 



— sin"" 1 



X 



/^ 3 + y 



2 



> 



x^ 



r 



The potential of the matter outside the cylinder *r~\2+/ m2 
.x. X y between the limits m=ju, and w=;L 



1 is the integral of 



m 



3 2 



a + b 



This gives 



/ 1 _JL\i a ~~^ 

u?] 2 a + b 
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Therefore the potential at the internal point x, y is 



c a — b cab . . , 

2 a + Z> u ' 2 » 



03 



V 



X% if 

a 2 W 1 



sin ] 



# 



V^ 3 + 2/ 3 



# 3 3/ 3 



Thus if the density at the point as, ^/ inside the elliptic cylinder — + — = 1 be 

a u 



e /1_1 



xy 



[a* by 



I, then the potential at an internal point is 



c a—h cab 

2a~+b Xy+ T 



S 



sin 



x 
a 



-> 



I 



v 



«# 3 ?/ 2 
a 3 Z> 3 



— sm 



-l 



# 



•\A 3 + f 



>, 



and the potential at an external point is 



a 



cab I a? -f Z> 3 + 2e \ c<z& J' . _, 



57 



2 p- ^8+ 



wher 



11/S 



2 



A i y 



sm 



-] 



^ 



\/^ + 2/ s 
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